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ABSTRACT 


During the month of July 1987 an acoustical experiment was 
conducted by the United States Naval Research Laboratory (NRL) 
in the East Greenland Sea Marginal Ice Zone (MIZ). Ambient 
noise "hot spots" or concentrated areas of relatively high 
noise levels were found along the ice edge using a towed 
array. Ambient noise levels were obtained on 27 and 28 July 
using AN/SSQ-57A and AN/SSQ-57XN5 calibrated sonobuoys. The 
temperature structure of the area was determined using XBT 
(ship) and AXBT (P3C aircraft) buoys placed inside and outside 
the ice edge. The ice edge was determined from coincident 
satellite photos, 90 GHz microwave imagery and P3 radar ice 
edge maps. Weather data (sea state and wind speed and direc- 
tion) were recorded on the ship. The data seem to indicate a 
correlation between the high ambient noise levels of the hot 
spots and the presence of a large topographically controlled 


mesoscale eddy located at the southeastern extent of the MIZ. 
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T INTRODUCTION 


The Arctic is a region of increasing interest to the 
United States Navy due its significant strategic value. 
Conducting Naval operations in this region requires an 
accurate knowledge of the environment and thus a primary goal 
of naval basic research in the Arctic is to develop a better 
understanding of Arctic processes [Curtin, 1988]. In July 
1987, an acoustic experiment was conducted by the United 
States Naval Research Laboratory (NRL) as part of operation 
RESOLUTE SUPPORT in the Greenland Sea. The purpose of the 
investigation, the first of two planned experiments with the 
Same objectives, was to identify and study local regions of 
high ambient noise levels or "hot spots" present in the 
marginal ice zone. An important aspect of the research was to 
examine the spectral characteristics of the hot spots by 
combining real time acoustic processing of the ambient noise 
with advanced remote sensing of the environment [NRL proposal, 
1987]. 

A marginal ice zone (MIZ) is a transition region from open 
ocean to pack ice where interaction between polar and temper- 
ate climate systems result in an edge of ice cover with strong 
horizontal and vertical gradients in the atmospheric and 
oceanic variables [Johannessen et al., 1983]. Hot spots, 


concentrated areas of relatively high ambient noise levels, 


have been shown to exist in the vicinity of the ice edge [Yang 
et al., 1987]. The hypothesis proposed to explain these local 
noise sources is that they are associated with mesoscale 
eddies embedded along the ice edge and are generated by two 
fundamental mechanisms: 1) collisions and crushing of ice 
floes in regions of eddy convergence and 2) breaking of ice 
floes in regions of eddy divergence [NRL proposal, 1987]. 
This thesis is the first comprehensive analysis of the 
data resulting from the first experiment. The second experi- 
ment was originally scheduled for Spring 1988, with the author 
to be intimately involved in the collection of data, but 
unfortunately the experiment had to be postponed until 1989. 
It will include more extensive data collection due in part to 
the lessons learned from this research and serve to validate 
the 1987 experimental results. Due to the time constraints 
involved in the completion of this thesis, the analysis 
presented herein is not meant to be final, but to serve as an 
in-depth preliminary look at the available data as well as 
data collection procedures. Hopefully, this will be beneficial 
to the planning, execution, and evaluation of the second 


experiment. 


CE- BACKGROUND 


A. PHYSICAL OCEANOGRAPHY 

The Greenland Sea is a very complex acoustic environment 
that until recent years had not received substantial attention 
from the scientific community. Thanks primarily to the 
Marginal Ice Zone Experiment (MIZEX) program, an Office of 
Naval Research (ONR) Arctic Sciences Program Accelerated 
Research Initiative (ARI) of five year duration [Curtin, 
1988], the knowledge base has increased significantly. This 
is readily apparent from the number of recent papers appearing 
in scholarly journals [see e.g., J. Geophys. Res. 92(C7), 
1987] and technical reports. 

In order to better understand the relationship between 
eddies in the marginal ice zone and the resultant ambient 
noise, it is necessary to first briefly examine the physical 
oceanography of the Greenland Sea. This cverview will focus 
on the primary features of the bathymetry and circulation, the 
water masses, the ocean frontal systems, and the marginal ice 
zone. The circulation in the Greenland Sea is dominated by a 
large cyclonic gyre formed by the southward movement of Polar 
Water (PW) from the East Greenland Current (EGC) north of Jan 
Mayen Island, the eastward flowing Jan Mayen Polar Current, 
and the northward and westward movement of part of the 


Atlantic Water (AW) in the north Greenland Sea [Johannessen, 


1986]. Aagaard [1970] has shown that the gyre is forced by the 
wind stress over the Greenland Sea but noted that thermohaline 
effects were also important. In addition, the circulation and 
distribution of the water masses is known to be strongly 
affected by the bottom topography [Johannessen, 1986; Gascard 
et al., 1988], e.g., the East Greenland continental shelf and 
a combination of fracture zones as seen in Figure 1. 

A representative temperature profile from outside the ice 
edge and its corresponding sound speed profile are shown in 
Figure 2. They are typical of cold Arctic water augmented with 
summer surface heating which results in a shallow sound 
channel axis. In general, the sound channel axis depth varies 
from O to about 100 m east of the ice edge. Under the ice and 
in the open sea in winter, the acoustic channel axis is at the 
surface producing the well-known "half-channel" (Figure 3). 
The major variation of axis depth is found along the Jan 
Mayen, Mohns, and Knipovich ridges as the axial depths rise 
from the southeast to northwest side of the basin. These 
locations are portrayed in Figure 4. 

Figure 5 shows the locations of the major fronts found in 
the Nordic Seas. These fronts are often anchored by bottom 
features such as continental shelf edges and midbasin ridges 
[Triantos and Hurdle, 1987]. The one of most significance to 
this experiment is the East Greenland Polar Front (EGPF). In 
summer it has a strong horizontal temperature gradient which 


on its cold water side rapidly brings the sound speed axis to 
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A map showing the general bathymetry 
100 m) and circulation in the Greenland Sea [from Paquette et 
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Figure 5. The major oceanographic fronts found in the Nordic 
seas [from Johannessen, 1986]. 


the surface. This front strongly affects acoustic propagation 
of surface noise sources by ducting energy from a surface 
channel under the ice downward into a 1000 m thick channel 
with an axial depth of about 550 m. The EGPF coincides with 
the high speed jet of the EGC which has recently been shown 
to be quite narrow (approximately 25-40 km), resemble a 
dynamically stable jet, and closely follow the continental 
slope [Gascard et al., 19881. The continental slope is 
understood to be the region between the 1000 and 2000 m 
isobaths. 

The marginal ice zone is a noisy, complex transition 
region which extends from the ice edge to the pack ice, about 
100 km into the ice covered region. In winter most of the ice 
in the Greenland Sea is comprised of thick first-year ice and 
multiyear ice. During the summer and fall, the majority is new 
or first-year ice. Gascard et al. [1988] characterize the MIZ 
below 77930'N in terms of two different limits between the 
pack ice and the open ocean. The first is the one at the sea- 
ice boundary which they term the outer limit and is commonly 
accepted to be the real ice edge. It is described as very 
diffuse in some places while being well defined by compacted 
small floes in others due to wind and swell resulting in a 
highly variable boundary. The second or inner limit deliniates 
the MIZ proper and the pack ice (ice concentration of 4/10 or 
greater) drifting in the EGC. It is more stable in location 


and more uniform than the outer edge but summer melting can 


10 


reduce its sharpness. The MIZ between these limits is char- 
acterized by low ice concentrations, highly variable width, 
typical floe size (diameter less than 1 km), and eddy fea- 
tures. Breakup of floes due to wave action and lateral melting 
contribute to this low concentration and floe size distri- 
bution up to 60 km into the ice [Gascard et al., 1988], but 


cannot account for the inner limit. They conclude that the 


inner limit must be related to the strong shear existing 
between the EGC central flow, where Arctic ice is moving 
southerly at speeds up to 1 kt, and the ice margin with its 
embedded mesoscale eddy circulation. During the month of July 
they observed strong eddy activity and a shift from northerly- 
northeasterly on-ice winds at the beginning of the month to 
southerly, off-ice prevailing winds late in the month. 
Mesoscale eddies found at the ice edge are often charac- 
terized by a pair of ice-tongue systems. The southerly one is 
due to off-ice entrainment (W to E) directly associated with 
cyclonic eddy rotation while one to the north is due to a 
divergence area created upstream [Gascard et al., 1988]. They 
reasoned that since the flow is basically quasi-irrotational, 
the velocity components wholly depend upon a velocity poten- 
tial function. Thus an area characterized by intense accel- 
eration such as is encountered near eddies should be counter- 
balanced by an area of deceleration and divergence to satisfy 
continuity. The eddies advect warm Atlantic Water westward 


beneath the ice (transporting heat into the MIZ) and pull ice 
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and Polar Water away from the ice pack into the warm AW. These 
features can often be observed by aircraft or remote sensing 
techniques. After ruling out the possibility of the EGC to 
self-generate ice-edge eddies, Gascard et al. [1987] found 
the interaction of westward propagating open ocean eddies from 
the West Spitsbergen coast with the EGC to be the prime eddy 
generating mechanism instead of fourth in the hierarchy of 


five sources given by Johannessen et al. [1987]. 


B. AMBIENT NOISE 

Ambient noise is the prevailing, sustained, unwanted 
background noise of the ocean. Ambient noise excludes self 
noise, electrical noise, and momentary sounds such as that of 
a nearby ship or biological organism (e.g., whale) [Urick, 
1986]. In essence, it is the remaining noise after all other 
noise sources have been eliminated. An excellent introduction 
to the topic of ambient noise, which includes material from 
a great number of ambient noise studies (through 1986), can 
be found in Urick's book, Ambient Noise In the Sea [1986]. 

Interest in the acoustic properties of the MIZ has 
increased manifold since the early 1970's when Diachok and 
Winokur focused their attention on the study of ambient noise 
along the Arctic ice/water boundary [Diachok and Winokur, 
1974]. Their significant finding was that ambient noise levels 
peaked at the ice/water boundary of the MIZ when compared with 


levels under the ice and in open water. They showed that this 
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peak was relatively sharp for a compact ice edge, becoming 
broader as the ice edge becomes more diffuse as illustrated 
in Figures 6a and 6b. Another feature evident in the two 
figures is the difference in the rates at which the ambient 
noise levels decrease with increasing distance from the ice 
edge. The levels under the ice fall off more rapidly than 
those in the open ocean. In addition to finding the ambient 
noise levels lower under the ice, they also found them to be 
intermittently spikey which was attributed to the collision 
and cracking of ice floes [Diachok, 1980]. Diachok concluded 
that sea-ice boundary-generated noise may be treated as a line 
source of sound at distances from the edge which are large 
compared to the magnitude of ice-edge irregularities [Diachok, 
1980]. 

Recent investigations by NRL scientists [Yang et al., 
1987] have shown that at the ice edge ambient noise is 
sometimes concentrated in small areas or "hot spots" along the 
ice edge. A feature which distinguishes hot spot noise from 
standard sea noise is a 12 dB/octave rolloff from 100 to 500 
Hz as opposed to a 6 dB/octave rolloff in the open sea. The 
following parameters are well known to affect ambient noise 
generation in the Arctic seas: sea ice kinematics (trans- 
lation, vorticity, divergence, normal deformation, shear 
deformation); sea state; currents; eddies; winds; air tempera- 
ture; and biologics [NRL, 1987]. The NRL scientists reported 


that they were unable to indentify the generating sources but 
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Figure 6. Variations of median ambient noise sound pressure 
spectrum levels with distance from the ice edge for 
frequencies of 100, 315, and 1000 Hz: (a) Compact ice edge 
with sea state 2; (b) Diffuse ice edge with sea state 1 [from 
Diachok and Winokur, 1973]. 
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the nominal separation (50 km) and spatial extent (5 km) of 
the noise hot spots were the same dimensions as some of the 
eddies found in the Greenland Sea by Johannessen et al. 
[1983]. What was needed was an experiment which combined real 
time acoustic processing with advanced remote sensing tech- 
niques to identify the primary contributors to ambient noise 
hot spots at low frequencies. This was the motivation behind 


the experiment reported herein [NRL proposal, 1987]. 
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TETS DATA COLLECTION AND ANALYSIS 


A. TYPES OF DATA 

To support the objectives of the field program data were 
collected from three entities: the research ship, aircraft, 
and satellite. Personnel on the ship recorded weather obser- 
vations, launched 47 expendable  bathythermograph (XBT) 
devices, and obtained navigational data and "hot spot" 
locations. A US Navy P3C (II) aircraft deployed 18 AN/SSQ-57A 
and AN/SSQ-57XN-5 ambient noise monitoring buoys and 25 AXBT 
buoys. An NRL RP3A aircraft gathered 90 GHz microwave imager 
remote sensing information. Satellite AVHRR (Advanced Very 
High Resolution Radiometer) data consisted of infrared (IR) 
and visual depictions of the ice edge and sea surface temper- 


ature (synoptic) gradients. 


B. DATA COLLECTION 
1. Ambient Noise Data 

Eight LOFAR sonobuoys (AN/SSQ-57A and AN/SSQ-57XN-5) 
were deployed on 27 July 1987 and ten on 28 July. LOFAR 
sonobuoys employ hydrophones which are omnidirectional in both 
the vertical and horizontal planes. The SSQ-57A is unique in 
that it is calibrated to allow accurate determination of 
ambient noise levels. The difference between a 57A and a 57XN- 


5 is that the latter has an operating depth of 1000 ft (305 m) 


16 


whereas the former has a maximum depth of 400 ft (122 m). The 
ambient noise data were collected on board P3C aircraft 
equipped with AN/ARR-72 receivers and a 28-track analog, FM 
wideband group II (WBII) tape recorder. The tape speed for 27 
July was 3-3/4 ips and 7-1/2 ips for 28 July. 

The ship, directed toward regions of high ice concen- 
tration by the aircraft and towing an acoustic array of 32 
hydrophones, moved parallel to the mean ice edge. If the ice 
edge acted as a continuous line source of constant source 
level, the noise would be distributed across many beams of the 
array. If, however, the ice edge noise was concentrated in 
localized "hot spots", then a peak would be evident in a small 
number of dominant beams [Yang et al., 1987]. The latter was 
found to be true and standard localization techniques were 
employed to fix the location of the hot spots. The array was 
used only to locate the hot spots. That is, the ambient noise 
analysis in this thesis does not include any array data. 

The buoys were deployed in pairs with one shallow (60 
ft) and one deep (300 ft or 1000 ft) in the vicinity of the 
hot spots. There were originally three data tapes, one for 27 
July and two for 28 July. Unfortunately, only the first of the 
two tapes for 28 July was retrieved after the experiment and 
thus some of the ambient noise data collected was not avail- 
able for analysis. Figure 7 shows the location of the ambient 
noise buoys in relation to the hot spots and radar ice edges 


for 27 and 28 July. 
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2. Bathythermograph Data 
Fifteen Airborne Expendable Bathythermograph (AXBT) 
buoys (AN/SSQ-36) were deployed on 27 July and ten on 28 July. 
Additionally, the ship deployed XBTs at midnight and every 
two hours between 0600 and 1800 on 24, 25 and 26 July. On 27 
July, a total of 12 were expended (additional buoys expended 
between 1200 and 1800). Eight were deployed on the 28th. 
Figure 8 shows the locations of the XBT and AXBT buoys in 
relation to the hot spots and radar ice edges. 
3. Weather Data 
Wind speed and direction and sea state were recorded 
every two hours from 24 July to 29 July by personnel on board 
the research vessel. 
4. Remote Sensing/Satellite Data 
An NRL RP3A equipped with a 90 GHz microwave radio- 
metry imager flew over the study area on 25 July and 28 July. 
AVHRR satellite data consisting of visual and (IR) imagery was 
collected from 24-27 July. Radar maps of the ice edge were 
made by ice observers on 27 and 28 July as shown in Figures 


7 and 8. 


C. DATA ANALYSIS 
1. Ambient Noise Data 
The 28-track tapes containing ambient noise and AXBT 
data were converted to 14-track tapes by Commander Task 


Group 168.2 in Norfolk, VA because a 28-track recorder was not 
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available at the Naval Postgraduate School (NPS). A Honeywell 
5600E Intermediate Band (IB) tape recorder was used to replay 
these tapes after performing initial calibration procedures. 
The Honeywell 5600E had only two available speeds: 7-1/2 ips 
and 15 ips. Since the 27 July tape was recorded at 3-3/4 ips, 
it had to be replayed at double speed. The 28 July tape was 
played back at the recorded speed of 7-1/2 ips. An inter- 
mediate band (IB) recording at 3-3/4 ips has a frequency 
bandwidth of only 1250 Hz while an IB recording at 7-1/2 has 
a bandwidth of 2500 Hz. The time code on the 27 July tape was 
not usable. This was true for the first copy of the 28 July 
tape as well. However, a 28 July tape with a usable time code 
and eight channels (as opposed to four) of ambient noise data 
arrived on 15 November 1988, too late to be of use in this 
thesis. 

The hard copy gram analysis records (flight grams) 
were used to determine what sonobuoy data was on each channel 
of the tape and at what times since the grams were hardwired 
to the recorder channels. Thus while different buoys could be 
on the same tape channel, a single gram was always associated 
with one particular channel of the tape. The grams were not 
as completely annotated as desired but after some effort a 
tape map was made delineating what was on each channel of the 
tape at any given time. 

Since a time code was not available, another means of 


determining times on the tapes was needed. The solution was 
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to select a channel with a known AXBT trace (easily identified 
on the gram), replay that channel from the beginning while 
listening for the start of the first AXBT and then to use it 
as a common reference point. The tape footage counter index 
at this tape position was assigned the AXBT time. It was then 
relatively easy to determine times and corresponding tape 
index values from a knowledge of the number of feet of tape 
travel and the tape speed. This procedure worked quite well. 
The tapes were played back many times before formally proces- 
sing them to select an interval free of RF interference, 
signal loss or anything else that would adversely affect 
processing. After the channels and their start/stop points had 
been identified using the times on the grams, the tapes were 
ready to be processed. Figure 9 shows the basic equipment con- 
figuration. 

A Spectral Dynamics Model SD375 Dynamic Analyzer II 
was used to analyze the sonobuoy data. This versatile instru- 
ment is a stand-alone, hard-vired, dual channel micro- 
processor-based FFT analyzer and signal processor that 
analyzes frequencies up to 100 kHz vith 400-line per channel 
resolution. An NEC APC IV personal computer (IBM/AT com- 
patible) equipped vith a National Instruments GP-IB interface 
card and softvare vas used to control all aspects of the data 
acquisition process. The control program vas vritten specific- 
ally for this thesis in Microsoft QuickBasic 4.0. The complete 


program listing is included as Appendix A. 
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After the descriptive information about the two 
channels, the frequency (frequencies) of interest and start 
time were entered, the tape was started (manually) and data 
acquisition begun. The sonobuoy data were linearly averaged 
for a period equal to 1 minute of real time data which equated 
to a 30 second processing interval for the 27 July tape and 
a 60 second interval for the 28 July tape. One minute averag- 
ing was chosen to average out the many temporary fluctuations 
and transients common to an ocean environment as well as to 
enhance identification of a persistent noise source or signal. 
The overall ambient noise power spectrum level was desired and 
not the spectra associated with transients. 

Hanning window amplitude weighting, chosen to minimize 
the effect of the window function on the analysis and a 
standard choice for ambient noise power spectral analysis, and 
maximum overlap processing (the default for any frequency span 
greater than 1 kHz) were used. The window function dramatic- 
ally reduced the sidelobe structure but had the unfortunate 
consequence of widening the resolution bandwidth (BW). For the 
Hanning window, the resolution bandwidth was computed from 
[Instruction Manual, 1981] 

BW - (Frequency Span/400)*1.5. 
Three different analyzer configurations were used as 
illustrated in Table I. Setups 1 and 2 were used to analyze 
the 28 July data. Two entirely separate analyses of the same 


section of the tape were done for the purpose of comparing the 
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Table I. SD 375 OPERATING PARAMETER SUMMARY 


Setup 1 Setup 2 Setup 3 
Weighting Hanning Hanning Hanning 
Frequency Span 0-2 kHz 0-400 Hz 0-4 kHz 
Averaging Linear Linear Linear 
Average Time 60 sec 60 sec 30 sec 
Input Level 2 Vrms 2 Vrms 1 Vrms 
Effective BW Zub. Hz Luiz 15 Hz 
Sample Rate 5.12 kHz 1.024 kHz 10.24 kHz 


ambient noise levels obtained using different frequency spans 
and thus resolution bandwidths. Setup 3 was used to analyze 
the 27 July data. Its frequency span and sampling rate were 
double that of Setup 1 to account for the frequency doubling 
resulting from playback at twice recorded speed. The averaging 
period was one-half of the Setup 1 value so that 1 minute of 
real time data would be represented by each sample. This 
produced basically identical processing of the two tapes even 
though they were recorded at different speeds. The sample rate 
was fixed at 2.56 * Frequency Span. Anti-aliasing filters were 
selected automatically by the SD375 on each channel for each 
frequency range with initial 120 dB/octave rolloff for 70 dB 


rejection of anti-aliasing terms [SD375 User's Manual, 1982]. 
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The computer checked the SD375 continuously near the 
end of the averaging interval until it detected that averaging 
Was complete. The SD375 then transmitted the precision data 
contents of its averager memory to the computer's random 
access memory (RAM) using ASCII byte coding which took about 
three seconds to complete. Immediately after receiving the 
data, the computer reset the averager memory and started the 
next sample. This ensured maximum utilization of the relative- 
ly small amount of useful data on the tape. Next, the channel 
À and channel B data were written to separate files on the 
computer's hard disk. The control phase of the computer 
program terminated when data sampling was stopped (manually) 
at the previously selected time. 

The remainder of the computer program handled file 
conversion. The ASCII files were converted to random access 
numeric files with one real number for each of the 400 bins. 
Finally, this number became the actual value of the averager 
memory for a single bin when multiplied by the input level 
selected (full scale value) and divided by 451 (the gain 
factor for a Hanning window [Instruction Manual, 1981]). 

The above procedure was repeated for each record in 
the input (ASCII) file resulting in a random access file 
containing 400 values for each sample conducted. To be useful 
this file needed to be subdivided into smaller files, one file 
for each sample. The program converted the random access file 


into ASCII files containing 400 records, one record for each 
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bin. In addition, an interactive file conversion subprogram 
created similar files for each frequency requested. These 
frequency files contained the values for a particular fre- 
quency from all of the samples. 

The statistical analysis of the data was accomplished 
with STATGRAPHICS [1988] statistical software. The files 
generated were imported directly into the program. Each file 
was converted into a variable array and all the arrays for a 
particular channel were grouped into a common file. The 
magnitude values were then transformed into corrected power 
values (from volts to dB re luPa?/Hz). 

Mathematically, this transformation was as follows: 

P = 20*LOG(M/lvolt) - 10*LOG(BW) + C(f) 
where P is the ambient noise level (power) in dB re 1uPa?/Hz, 
M is the magnitude of each bin in volts, BW is the effective 
bandwidth of each bin and C(f), a function of frequency only, 
is the correction for each bin frequency which relates dB re 
1V2 to dB re luPa?. The standard sonobuoy calibration curve, 
with a published uncertainty of +2 dB, used to determine the 
values of C(f) was provided by NRL [Gruber, 1988] but is clas- 
sified and therefore will not be discussed in this thesis. 
There is an additional correction term, not shown in the above 
equation, involved which accounts for the receiver calibra- 
tion. This term was omitted because it was unable to be 
determined from the data tapes. That is, known calibration 


tones were not available for each channel. This has an obvious 
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effect on the accuracy of the results. This correction is 
assumed to be constant with frequency (commonly done since 
only one frequency of known amplitude is usually used) [NRSC, 
1988]. If the reference value of 1 volt was actually 1.1 
volts, then this would amount to a correction of -0.8 dB. That 
is, the true value would be 0.8 dB lower than the values 
reported herein. 

The individual spectra of each sample were averaged 
to create a mean spectrum for each buoy during the time frame 
investigated. Averaging was performed on the magnitude data 
to preclude errors due to nonlinearities from entering into 
the data analysis. These mean magnitude data were then trans- 
formed into corrected mean power spectra in the method 
described above and plotted. The mean power spectra for each 
buoy is included in Appendix B. 

À more descriptive and statistically sound analysis 
was performed by determining the 50th percentile value 
(median) of all the values in the set of samples for each fre- 
quency of interest. The spread in the ambient noise measure- 
ments at each frequency was determined from the difference 
between the 90th and 10th percentiles. By working with the 
median vice average value, dramatic shifts in the magnitudes 
of a few samples would not have undue influence on the value 


sought to be the most representative of the data set. 
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2. Oceanographic/Environmental Data 
a. Bathythermograph Data 
As an aid to making comparisons between the 
ambient noise levels of different buoys, a knowledge of the 
local environmental differences in the conditions pertaining 
to each buoy is very important. For this reason temperature 
contour plots (Figures 10, 11, and 12) were drawn for depths 
of 16, 33, and 60 ft (5, 10, and 18 m), respectively. The 
60 ft depth corresponds to the depth of the shallow buoys. The 
deep buoys on 27 July were placed at 1000 ft while on 28 July 
they were located only 300 ft below the surface. All XBTs were 
treated synoptically, i.e., as if they were expended at the 
same time. 
b. Weather 
The wind speed and direction and sea state 
recorded every two hours were plotted (Figures 13 and 14) to 
determine not only the conditions which existed during the 
time frame of sonobuoy data acquisition but also to examine 


how these conditions changed over time. Table II is a summary 


Table II. WEATHER DATA SUMMARY 
SSS n. == aamoaamMoooooooo 


24 July 25 July 26 July 27 July 28 July 


Sea State 0-1 1-4 I ls OE 
Wind Speed (kts) 4-8 3-20 6-10 5-15 0-4 
Wind Direction (27) 280 220 240 275 230 
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the range of conditions on 24, 25, 26, 27 and 28 July. 
These parameters have been known to fluctuate dramatically 
over relatively short periods of time. 
c. Remote Sensing and Satellite Data 

Radar maps of the ice edge were made by P3 person- 
nel on 27 and 28 July. All of the 27 and 28 July ice edges 
appearing in the figures are based on these radar maps. 
Microwave radiometer imagery at 90 GHz was obtained from a 
high altitude ice-edge overflight by NRL's RP3A on 25 July. 
The effort on 28 July did not produce useable results. The 
satellite information was processed by NRL and forwarded to 
NPS. This consisted of: xerox copies of seven NY satellite 
photos (5 visual and 2 IR for the period 24-27 July); four 
transparencies made from the visual satellite photos which 
were enlarged and cropped to show the operational area (one 
for each day); and two color transparencies made from the IR 
photographs (24 and 25 July). The xerox copies can be found 
in Appendix C. The most comprehensive and useful information 
available was for 25 July. No satellite information was 
obtainable for 28 July. The 25 July color IR transparency was 
used to trace out the ice edge and surface temperature 
contours (Figure 15). Very much enlarged, the important 
portion of Figure 15 is shown in Figure 16. None of the satel- 
lite data items contained a latitude/longitude grid. This made 
it difficult to derive the exact position of the ice edge in 


relation to the sonobuoy and XBT positions from the satellite 
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Figure 15. Sea surface temperature (SST) gradients extracted 
from a satellite color IR photo for 25 July 1987. The hatched 
line indicates the ice edge. The warm water is about 6°C. The 
satellite photo did not contain geographic coordinates. 
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Figure 16. Enlargement of Figure 15 to isolate the 
operational area (indicated by the box) of the experiment. 
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photos alone. However, a fairly accurate estimate was made by 
relating known geographic features from a map to those on the 
small scale photographs and estimating the coordinates of the 
ice edge. The 90 GHz imager pictures were helpful also since 
the start and finish latitude and longitude of each leg were 
known. Figure 17 shows a black and white version of the 
original color print provided by NRL reduced in size. The 
latitude and longitude are approximations since the six legs 
were not flown at exactly the same height and orientation and 


thus are also at slightly different scales. 
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IV. RESULTS 


A. AMBIENT NOISE DATA 
1. 27 July 1987 

Of the eight channels used for ambient noise data 
collection, only four could be analyzed since the hard copy 
gram (each gram has four channels recorded on it) for the 
other four channels was not retained after the experiment. The 
four channels which did have a corresponding gram contained 
ambient noise data for five buoys. All buoys are referred to 


by their transmission channel number. Table III summarizes 


Table III. MEDIAN AMBIENT NOISE LEVELS 27 JULY 87 





Frequency Buoy 4 Buoy 6 Buoy 12 Buoy 16 Buoy 23 


30 Hz 78 76 71 73 70 
40 Hz 80 19 13 19 12 
50 Hz 82 19 74 TS 71 
100 Hz 81 13 76 80 73 
200 Hz 19 76 76 11 73 
315 RZ 78 74 72 74 70 
400 Hz 74 72 70 72 68 
500 Hz 15 70 69 70 66 
600 Hz 71 67 67 69 63 
700 Hz 69 66 67 67 62 
800 Hz 68 64 63 66 61 
900 Hz 67 62 63 66 60 
1000 Hz 68 63 62 66 59 
11009H2 66 61 61 65 58 
1200 Hz 65 60 60 64 56 


(Note: ambient noise levels are in dB re luPa?/Hz) 
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the median ambient noise levels recorded on 27 July for 
buoys 4, 6, 12, 16 and 23. Buoys 4 and 6 were set deep (1000 
ft) while 12, 16 and 23 were set shallow (60 ft). The number 
of samples (minutes) of analyzed data used to determine the 
median noise level was as follows: 17 for buoy 4; 14 for buoy 
6; 15 for buoy 12; 11 for buoy 16; and 32 for-buoy 23. The 
spread in the median ambient noise levels of Table Ill, 
determined from the 10th and 90th percentiles, is shown in 
Table IV. It is apparent that there is temporal variability 
in the ambient noise levels over the relatively short obser- 
vation period. 

Figure 18 illustrates how the spectral amplitudes of 


each of these buoys compare. The frequency range extends only 


Table IV. SPREAD IN AMBIENT NOISE LEVELS 27 JULY 87 
 — Im nuA— m | 


Frequency Buoy 4 Buoy 6 Buoy 12 Buoy 16 Buoy 23 


30 Hz Sr? 4.0 2.8 2.8 Dd 

40 Hz 2.6 2 1:9 4.5 T 

50 Hz 33 1 7.4 338 51 
100 Hz 254 2.5 245 179 1.5 
200 Hz 1 1.4 3.0 27 2:2 
%15. HZ 253 1.4 4.0 351 3.7 
400 Hz 3.0 2.3 251 3.8 21 
200 Hz JT 6.0 2:9 4.6 20 
600 Hz 2.8 1.6 30 5: 2.1 
700 Hz 4.3 1.8 2.6 57 21 
800 Hz 3.6 251 3.0 23 1.8 
900 Hz 6.8 1.6 7.8 2.6 2.6 
1000 Hz 9.5 0.9 1.0 5.8 2.8 
1100 Hz 7.0 1.5 1.8 1.4 2.4 
1200 Hz 4.9 1-6 1.4 2 2.4 


(Note: spread in ambient noise levels is in dB re luPa?/Hz) 
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to 1200 Hz due to the smaller bandwidth of the recording as 
discussed earlier. Appendix B contains a separate plot of the 
median noise spectrum levels for each buoy and the associated 
10th and 90th percentile values. Included in Appendix B are 
separate plots of the mean ambient noise levels. From Figure 
7, the loudest buoys (4 and 16) were located at the ice edge 
and close to the hot spot as determined by triangulation from 
the research vessel. Buoy 6, also exibiting significant noise 
levels, was near the previous day's hot spot and surrounded 
on three sides by ice at a distance of about 8 km. Buoy 23, 
the quietest buoy, was approximately 2 km south of the ice 
edge. Buoy 12, on the other hand, had levels in between those 
of 16 and 23 and was located just inside the ice edge to the 
east which forms a sort of penninsula feature. 
2. 28 July 1987 

The 14-track tape copy analyzed for this thesis 
included only four of the original eight ambient noise data 
channels. The second copy of this day's tape which arrived on 
15 November 1988 contained all eight sonobuoy channels as well 
as a usable time code, but was too late to be of use as stated 
before. For the 35 minute interval of the tape studied, each 
of the four channels contained ambient noise data for a 
separate buoy. The buoys analyzed were 18, 20, 28, and 29. 
These buoys account for two of the four pairs that were 
deployed. The first pair was comprised of 18 (shallow, 60 ft) 


and 28 (deep, 300 ft) which were directly on top of the hot 
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spot about 8 km west of the ice edge. The second pair, 20 
(shallow) and 29 (deep), were placed to the northwest of the 
first pair and approximately 8 km from the ice edge. 

Two analyses were done on this data as noted pre- 
viously: one with an effective resolution bandwidth of 1.5 Hz 
and one with 7.5 Hz. The 7.5 Hz resolution bandwidth, cor- 
responding to a frequency span of 0-2000 Hz, was the basis for 
Setup 1 as depicted earlier in Table I. The 1.5 Hz high 
resolution bandwidth, corresponding to a frequency range of 
0-400 Hz, was the basis for Setup 2. To distinguish between 
the different bandwidths for the same buoy in the discussion 
and tables which ensue, a buoy number followed by (A) indi- 
cates Setup 1 and when followed by (C) indicates Setup 2. 
Table V summarizes the median ambient noise levels for both 
setups. Figure 19 shows graphically how the median ambient 
noise levels of all the buoys determined by using Setup 1 
compare. The median ambient noise levels resulting from use 
of Setup 2 are shown in Figure 20a. To illustrate how the 
results from the two different setups compare, the results 
from Setup 1 in the range O to 400 Hz are shown as Figure 20b 
on the same page. It is readily apparent that they compare 
very favorably. The spread of each of the ambient noise levels 
listed in the previous table are given in Table VI. The 
temporal variability of the ambient noise levels noted earlier 
for 27 July is evident on this day as well but to a slightly 


higher magnitude. 
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Figure 20. Comparison of ambient noise levels (power spectra) 
from 0-400 Hz for buoys deployed on 28 July 1987 obtained with 
different resolution bandwidths. (a) Setup 1 with an effective 
resolution of 7.5 Hz. (b) Setup 2 with an effective resolution 
Of 1 59H27. 
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Table V. MEDIAN AMBIENT NOISE LEVELS 28 JULY 87 
A r sc "—————————— ЫЫЫ 


Frequency Buoy 18 Buoy 20 Buoy 28 Buoy 29 
(A) (C) (A) (C) (A) (C) (A) (C) 
10 Hz 95 97 104 104 105 105 101 104 
20 Hz 97 97 100 102 101 101 98 98 
40 Hz 95 93 98 97 98 98 95 95 
50 Hz 95 93 97 97 97 97 95 95 
100 Hz 88 88 93 92 93 93 93 92 
200 Hz 83 83 89 87 89 90 88 89 
515 Hz 78 78 85 82 85 84 84 85 
400 Hz 79 77 84 84 84 83 82 83 
500 Hz 76 80 83 81 
600 Hz 75 78 80 79 
700 Hz 73 76 719 77 
800 Hz 71 74 771 76 
900 Hz 70 73 76 74 
1000 Hz 73 713 74 74 
1100 Hz 68 71 74 12 
1200 Hz 69 69 73 71 
1300 Hz 66 68 71 70 
1400 Hz 65 69 70 70 
1500 Hz 64 66 69 68 
1600 Hz 63 65 68 67 
1700 Hz 62 64 68 67 
1800 Hz 62 64 67 67 
1900 Hz 61 63 66 66 
2000 Hz 61 62 65 65 


Note: Ambient noise levels are in dB re 1uPa?/Hz 
(A): Setup 1 [0-2 kHz] (C): Setup 2 [0-400 Hz] 
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Table VI. 
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B. OCEANOGRAPHIC/ENVIRONMENTAL DATA 

Figure 21 shovs the detailed bathymetry of the study 
region from 709N-749N, 09V-259V. In this figure, a significant 
easterly excursion of the 2000 m isobath, and thus the 
continental slope, occurs at about 72°30'N, 15°30'W associated 
with a long and narrow ridge which rises to a depth shallower 
than 1000 m. This significant bathymetric change diverts the 
EGC causing a portion of it to be shunted eastward and form 
the Jan Mayen Polar Current (JMPC) as shown earlier in Figure 
1. A more detailed view of this ridge is shown in Figure 22 
from which much of it is observed to lie within the oper- 
ational area (indicated by the grid) of the experiment. There 
is a strong suggestion that bathymetric steering of the EGC 
and ice edge occurs because of the presence of the ridge. 

The area under investigation lies at the southeastern 
corner of the local marginal ice zone as shown in Figure 23. 
In this figure, which is derived from Figure 15, the shaded 
regions indicate depths shallower than 1200 m and are esti- 
mates of the locations of the two ridges identified in the 
previous figure. The circled numbers represent two outflow 
regions of polar water. The first is believed to be due to 
that part of the EGC diverted eastward by the topographic high 
described above and also that part steered eastward by the 
smaller ridge (71°54'N, 13°45'W). The main contributor to this 
outflow appears to be the smaller ridge. The flow is initially 


eastward north of the main ridge and then southeastward as it 
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passes over the ridge and comes under the influence of the 
smaller ridge. The second region is possibly due to the south- 
eastwardly flow south of the smaller ridge and the Maro Bank 
region of the Jan Mayen Fracture Zone. This flow does not 
contribute to the JMPC. The patch of warm water in the lower 
right of the figure appears to lie over the Maro Bank region. 
The proposed flow pattern in the first region suggests a 
cyclonic circulation pattern. Figure 24 summarizes the 
hypothesised flow pattern. No current measurements were made 
to verify this flow but, the 5, 10, and 18 m temperature 
contour plots (Figures 10-12) support this hypothesis. 

The lack of geographical coordinates on the satellite 
photos complicates the analysis. In addition, the sea-ice 
boundary shows up very faintly, if at all, on the satellite 
pictures due to the resolution at which they were processed. 
From the four day sequence of satellite photos, the ice pack 
and MIZ, already quite diffuse, appear to have become more 
diffuse in the north-south direction as they moved slowly 
south with the currents but slightly more compact in the east- 
west direction. This is probably a result of the winds which 
blew from out of the W-SW and ranged in speed from 3 to 20 kts 
during the period 24-27 July before subsiding on 28 July to 
maximum of 4 kts. The maximum sea state went from 4 on 25 July 


277 on 28 July. 
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VIE. DISCUSSION OF RESULTS 


A. THE PRESENCE OF AN MIZ EDDY 

Although the remote sensing information does not explic- 
itly reveal the presence of an eddy, it does not discount its 
possibility either. The temperature contour plots (Figures 
10-12) in conjunction with the detailed bathymetry seem to 
reveal the presence of a large, topographically triggered 
mesoscale eddy centered at about 72°11'N, 012°50'W which is 
in the middle of the cyclonic circulation indicated earlier 
in Figure 24. It is estimated to be elliptical in shape, with 
major and minor axes of 50 and 30 km respectively, and 
cyclonic in circulation. It is oriented with its major axis 
running northwest to southeast. The position of the center of 
this proposed eddy with respect to the bathymetry, hot spots, 
and ice edge is shown in Figure 25 as located above a topo- 
graphic depression (with a depth greater than 2200 m) that 
lies between the two ridges. That portion of the ice edge 
observed over the four day period of this experiment is 
thought to have moved through the eddy from west to east while 


being steadily influenced by west-southwest winds. 
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B. 27 JULY 1987 

This day's hot spot was located in a region of ice 
divergence (near the beginning of cold water tongue, see e.g. 
Figure 10). This divergence is probably due to a combination 
of the presence of the eddy and the diversion of some of the 
EGC flow which aided the formation of the eddy. The proposed 
eddy location is east of this hot spot position. The small ice 
tongue apparent in Figure 7 is indicative of divergence and 
most likely the result of entrainment of the ice by the 
cyclonic circulation of the eddy. On-ice winds, averaging 9 
kts from 275°T, appear to produce an area of local convergence 
on the west side of the ice edge next to this hot spot. As can 
be seen from Figure 10, the strongest temperature gradients 
are to the immediate west and south of the hot spot. Thus the 
ambient noise levels observed are believed to be due to the 
combination of the breaking of ice associated with divergence 
and the ice collisions due to wind-induced convergence, with 
divergence being the dominant factor. Swell data would be 
helpful in this analysis of generating mechanisms but was not 
recorded during the experiment. 

Taking a look at the individual ambient noise levels, buoy 
4, located at the hot spot at a depth of 1000 ft (305 m), had 
the highest amplitdes at all frequencies. Buoy 16, also at the 
hot spot but at a depth of 60 ft (18 m), was the second 
loudest at almost all frequencies. Because of the config- 


uration of the ice and the assumed water flow, it is not 


57 


surprising that the noise amplitudes of these two buoys were 
louder than the other three. This lends creedance to the 
evaluation of the location as a hot spot. It is interesting, 
however, that the deep buoy was louder (1-3 dB) than the 
shallow one. This result is true for the other buoy pair (6 
deep and 23 shallow) as well. This latter buoy pair mwas 
separated by about 3 km and thus the two buoys were not as 
close together as the other pair. Buoy 23, with the lowest 
levels of all buoys on this day, was closer to the ice edge 
than buoy 6. The one remaining buoy, number 12, was located 
just inside (2 km) the ice edge almost 50 km northeast of 
buoys 6 and 23. Its levels were in between those of the other 
two shallow buoys. The observed levels were of the same 
magnitude as those of Johannessen et al. [1988] for 19 June 


84. 


C. 28 JULY 1987 

The hot spot identified on this day was approximately 2 
degrees of longitude (70 km) to the east of the 27 July hot 
spot. It is believed to be in a region of eddy convergence as 
shown in Figure 25. Also, from Figure 7, the large penninsula 
feature developed curvature on its west side as it moved to 
the east which supports the presence of the proposed eddy. 
This day was quieter than the previous day with respect to sea 
state and winds. However, the ambient noise levels were still 


very high for all frequencies. This same interesting result 
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was found by Johannessen et al. [1988] for 30 April 85. In the 
absense of strong winds or high sea state, the ice interac- 
tions are predominately caused by horizontal shear velocity 
arising from the current and eddy circulation [Johannessen et 
al., 1988]. The slight winds present were on-ice at the ice 
edge closest to the hot spot. This suggests that ambient noise 
levels in regions of eddy convergence may be of higher 
amplitudes than those in the regions of divergence due to the 
larger area and higher ice concentration of the covergence 
region. 

For buoy pair 18 (60 ft or 18 m) and 28 (300 ft or 91 m) 
positioned at the hot spot, the deep buoy is again louder at 
all frequencies. This holds true for the second pair as well. 
Buoy 20 (300 ft) is slightly louder than buoy 29 (60 ft) at 
all frequencies except 400 Hz where a signal of unknown origin 
is believed to have affected the data. The shallow buoys were 
expected to be louder than the deep buoys since they were 
located in the surface sound channel where the ice interac- 
tions were occuring. Distant heavy shipping may account for 
the fact that the deep buoys were louder but this is not known 
to actually be the case. There was no indication of the bouys 


malfunctioning. 
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VIT. CONCLUSIONS AND RECOMMENDATIONS 


In July 1987, NRL conducted an acoustic experiment in the 
Greenland Sea marginal ice zone (MIZ). Several ambient noise 
hot spots or concentrated areas of relatively high noise 
levels were found at the ice edge using a towed array. Ambient 
noise levels were obtained on 27 and 28 July with AN/SSQ-57A 
and AN/SSQ-57XN5 calibrated  sonobuoys. The temperature 
structure of the area was determined using XBT (ship) and AXBT 
(P3C aircraft) buoys placed inside and outside the ice edge. 
The ice edge was determined from coincident satellite photos, 
90 GHz microwave imager data, and P3C radar ice edge maps. 
Wind speed and direction and sea state were recorded on the 
ship. 

The study area was sampled well enough to indicate that 
a variable ambient noise field of significantly high levels 
existed in the region of the hot spots. In addition, the 
temperature contour data, bathymetry, and ice edge features 
and locations seemed to reveal the presence of a large 
mesoscale eddy embedded in the ice edge which is topo- 
graphically controlled. There appears to be a definite 
relationship between the observed hot spots and the presence 
of this eddy. The 27 July hot spot is believed to be as- 
sociated with a predominately divergent region (there was some 


local wind-induced convergence) while the 28 July hot spot is 
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probably associated with a region of convergence. The ambient 
noise levels on 28 July, the quietest day with respect to sea 
state and wind speed, were higher than those of 27 July 
suggesting that eddy convergence regions are louder than eddy 
divergence regions due to the larger area and higher ice con- 
centrations found in the former region. In the absense of high 
winds or sea state, the ice interactions are predominately 
caused by horizontal shear velocities arising from the current 
and eddy circulation. For every pair of colocated buoys 
analyzed, the deep buoy had slightly higher amplitudes. This 
was the opposite of what was expected. No explanation other 
than the possibility of distant shipping can be given to 
support this result. A larger number of buoys analyzed over 
much longer periods of time are needed to investigate this 
occurance. 

There are a number of recommendations, based on the 
analysis of this experiment, that should prove valuable to the 
follow-on experiment scheduled for next year. More oceano- 
graphic/environmental data should be taken such as CTD, 
current measurements, swell height and direction, detailed 
visual observations of the weather and ice, air temperature 
and barometric pressure. The bathymetry of the operation area 
should be a prime consideration in the placement of sonobuoys 
and XBT's since the visual observations will not always give 


accurate locations of the fronts and eddies. It would also be 
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helpful to place a buoy in a known "quiet spot" for comparison 
purposes. 

The procedures for collection of ambient noise data 
onboard the aircraft should be standardized for the experiment 
to facilitate analysis of the data. In particular, calibration 
tones of known amplitudes should be recorded on each channel 
of the tape. The flight grams need to be annotated to the 
highest degree of detail possible. A separate log (possibly 
just a pocket tape recorder) detailing buoy information, 
channel changes, etc. should be maintained. Most importantly, 
all of the data collected needs to be accounted for after the 


flight (correct number of tapes, etc.). 
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APPENDIX A 


QUICKBASIC 4.0 COMPUTER—CONTROL PROGRAM 


"PROGRAM: SD375.BAS 

"DATE: 10/19/88 

"AUTHOR: LT KRIS BIGGS, USN 

"PURPOSE: CONTROL SD375 DATA ACQUISITION AND FILE CONVERSION USING 
i NATIONAL INSTRUMENTS GP-IB INTERFACE CARD (IEEE-488) AND NEC 
i APC IV PERSONAL COMPUTER (IBM/AT COMPATIBLE). 

COMMON SHARED IBSTA%, IBERR%, IBCNT% 

DECLARE SUB InitialSetup () 

DECLARE SUB GetPrecisionData () 

DECLARE SUB PrintFiles () 

DECLARE SUB AskIfPrintoutWanted () 

DECLARE SUB CheckStatus () 

DECLARE SUB ResetClock () 

DECLARE SUB ConvertRawData () 

DECLARE FUNCTION MakeReal! (B1%, B2%, B4%) 

DEFINT A-Z ‘all untyped var are integers 

OPTION BASE 1 

TYPE RecordTypel 


DataDate AS STRING * 10 ‘date of data 
ChannelID AS STRING * 3 MA 2B etc 
TimeO AS STRING * 8 "SAMPLE TIME 
Freqi AS SINGLE '40 HZ 
Freq2 AS SINGLE '100 HZ 
Freq3 AS SINGLE "945 HZ 
Freq4 AS SINGLE "1000 HZ 
Freq5 AS SINGLE '1500 HZ 

END TYPE 

TYPE RecordType2 
Freq AS INTEGER 'center freq of bin 
Voltage AS SINGLE 'magnitude of voltage in bin 
Level AS SINGLE 'magnitude in dB 

END TYPE 


DIM ChA AS RecordTypel 

DIM ChB AS RecordTypel 

DIM Spectrum AS RecordType2 

DIM PrecisonData AS STRING * 6402 

DIM RawData AS STRING * 3200 

DIM Averaging%(1 TO 2) 

DIM Volt(1 TO 400) 

DIM Level(1 TO 400) 

CLS : SCREEN 0, 1: COLOR 14, 0: WIDTH 80, 25 

'this converts the precision data file created by a previous run PRINT 
"Do you want to do file conversion now [Y/N]?": AS = INPUT$(1) 
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IF UCASES(AS) = "Y" THEN 
CALL ConvertRawData 
CLS 
END IF 
'Data collection? 
PRINT "Do you want to collect data now (Y/N)?": AS - INPUTS (1) 
IF UCASES(AS) = "N" THEN GOTO EndProgram 
Filel$ = "Chana": File2$ = "ChanB" 
CALL IBFIND("SD375", SD375%): CLS 'open device 
CALL InitialSetup: CLS : COLOR 14, O 
PRINT TAB(19); "SD 375 is correctly configured and ready to begin.": 
PRINT 
PRINT TAB(19); "The front panel is not locked out.": PRINT 
COLOR 9, 0: LOCATE 10, 20: PRINT "Do"; : COLOR 30, 0: PRINT " NOT "; 
COLOR 9, 0: PRINT "change any of the front panel settings!" 
COLOR 14, 0: LOCATE 18, 20 
PRINT "Press [C] to continue or [A] to abort program": AS = INPUTS (1) 
WHILE NOT (UCASES(AS) = "A" OR UCASES(AS) = "C") 
LOCATE 18, 20 
PRINT "Press [C] to continue or [A] to abort program": AS = INPUTS(1) 
WEND 
IF UCASES(AS) = "A" THEN GOTO EndProgram ELSE CLS 
LOCATE 15, 25: PRINT "Continuing...." 
'open data files 
OPEN "ChanA.DAT" FOR RANDOM AS $2 LEN LEN(ChA) 
OPEN "ChanB.DAT" FOR RANDOM AS %3 LEN LEN (ChB) 
OPEN "LPT1:" FOR OUTPUT AS #4 'printer 
OPEN "MagA.DAT" FOR APPEND AS $5 LEN = 3200 "rav data chA 
OPEN "MagB.DAT" FOR APPEND AS #6 LEN 3200 'raw data chB 
RecordNumber2 = LOF(2) V LEN(ChA) 
RecordNumber3 = LOF(3) \ LEN(ChB) 
CLS ; COLOR 9, 0: LOCATE 3, 1 
PRINT "ChanA.DAT will contain the data for channel A": PRINT 
PRINT "ChanA.DAT currently contains "; NunberOfRecords2; " records." 
COLOR 15, 0: LOCATE 7, 1 
PRINT "ChanB.DAT will contain the data for channel B": PRINT 
PRINT "ChanB.DAT currently contains "; NumberOfRecords3; " records." 
LOCATE 11, 5: PRINT : COLOR 14, 0 
PRINT "The existing records in these files will be retained. Even so," 
PRINT "These files should be backed up to a floppy disk." 
PRINT : COLOR 25, 0: LOCATE 20, 18 


PRINT "---> Press any key to continue <---": COLOR 14, 0 
DO: LOOP UNTIL INKEYS <> "" 'wait for user to press a key 
DateToday$ = DATES 'saves actual date 
' k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k Kk 
DO 
CES 


LOCATE 8, 15 

INPUT "Date of Data tape (mm-dd-yyyy)? ", DataDate$ 

PRINT : PRINT "Start time of data sampling (hh:mm:ss)? "', StartTime$ 
DO: LOCATE 12, 15: INPUT "The Hour is (00 - 24) "; Hour$ 

LOOP UNTIL VAL(Hour$) >= O AND VAL(Hour$) <= 24 
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PRINT 

DO: LOCATE 12, 15: INPUT " and the minutes are (00 - 59) ".Minutes$ 
LOOP UNTIL VAL(Minutes$) >= 0 AND VAL(Minutes$) <= 59 

PRINT 


DO:LOCATE 12,15:INPUT" and the seconds are (00 - 59) ". Seconds$ 
LOOP UNTIL VAL(Seconds$) >= 0 AND VAL(SecondsS) <= 59 
StartTime$ = Hour$ + ":" + Minutes$ + ":" + Seconds$ 


LOCATE 12, 15: PRINT "The time entered is "; StartTimeS; SPACES (20) 
LOCATE 15, 15: PRINT "Are these correct (y/n)? ": AS - INPUTS(1) 
LOOP UNTIL UCASES(AS) - "Y": CLS 


DO 

LOCATE 10, 15 

INPUT "Enter the channel ID (i.e. 1-A) for CH A: ", ChA.ChannelID 
PRINT : LOCATE 12, 15 

INPUT "Enter the channel ID (i.e. 1-A) for CH B: ", ChB.ChannelID 


PRINT : LOCATE 15, 15 

PRINT "Are these correct (y/n)? ": AS - INPUTS (1) 
LOOP UNTIL UCASES(AS) = "Y": CLS 
ON TIMER(1) GOSUB DISPLAYTIME 
TIMER ON 


—“ mam l ea нән = —— — —— — =a a 
PRINT "Press [S] to start sampling, [E] to end.": AS = INPUTS(1) 
TIMES = StartTime$ 
DO UNTIL UCASES(AS) » "S" OR UCASES(AS) = "E" 

LOCATE 14, 20: PRINT "Incorrect response. Please try again.": BEEP 

AS = INPUTS(1) 

LOOP: CLS 

LOCATE 15, 30: PRINT "Conducting first sample": COLOR 1, O 
LOCATE 25, 50: PRINT "Press [F1] to quit sampling.": COLOR 14, 0 
IF (UCASES(AS) = "E") THEN GOTO EndProgram 

Averaging%(1) = 1 

ON KEY(1) GOSUB StopSampling 


KEY(1) ON '[F1] stops sampling 
CALL IBSAD(SD375%, 115) 

CALL IBWRT(SD375%, "MPB(A)/" + CHRS(13)) "MAKE PUSH BUTTONS ASCII 
CALL IBSAD(SD375%, 111) "PUSH BUTTON GROUP 15 
CALL IBWRT(SD375%, "39474846FF") "Mİ, STOP, ERASE, START 


TIMEŞ - StartTime$ 
Interval! = 59! 
T1! = TIMER 
DO: LOOP UNTIL TIMER - T1! > Interval! 'waits 59 sec to check status 
DO WHILE Averaging%(1) = 1 
CALL CheckStatus 
LOCATE 1, 1: PRINT "Checking Status...Average in progress..." 


LOOP 'keeps checking status until not 
averaging 
LOCATE 1, 1: PRINT "Averaging complete"; SPACES(50) 
CALL GetPrecisionData 'data for first sample 
DO 

T1! - TIMER 


DO: LOOP UNTIL TIMER - Ti! » Interval! - 1! 'waits 59 sec 
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LOCATE 1, 1: PRINT "Averaging complete"; SPACES(50) 
DO WHILE Averaging%(1) > 0 
CALL CheckStatus 
LOCATE 1, 1: PRINT "Average in progress..." 
LOOP "keeps checking status until not averaging 
LOCATE 1, 1: PRINT SPACES(70) 
CALL GetPrecisionData 
LOOP ‘loops until [F1] is pressed 
StopSampling: 
PrintOut: 
CALL ResetClock "reset system clock 
" CALL ConvertRavData 
EndProgram: 
" CALL AskIfPrintoutWanted 
CLS : PRINT "Program terminated." 
END 


t k x x k k k k k k k k k k k k k k k k k k K k k K k k K k K k k KX k k K k k k k k k k k k k k k k k k k k k k k k k k k k X k k 
t k k x k k k k k k k k k k k k k K END OF PROGRAM Kk k k k k k k k k k k k k k k k k k k k k k k k k k k k k k K 
t k k k k K k k k k k k k k k k k k k k k k k k k k k k k K k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k K k Kx 


DISPLAYTIME: 

LOCATE 2, 50: PRINT TIMES 

LOCATE 3, 50: PRINT SampleTimeS; " = last sample" 
RETURN 

SUB InitialSetup 'takes about 1 sec to execute 
SHARED SD375% 

CLS 

PRINT TAB(30); "FREQUENCY SELECTION: ": PRINT 
PRINT TAB(30); " (1) 40 Hz": PRINT TAB(30); " (2) 100 Hz" 
PRINT TAB(30); " (3) 315 Hz": PRINT TAB(30); " (4) 1000 Hz" 
PRINT TAB(30); " (5) 1500 Hz": PRINT TAB(30); " (6) All of then" 
COLOR 9, O 

LOCATE 12, 26: PRINT "Enter desired frequency (1-6)"; : F% = 


VAL(INPUTS (1)) 
WHILE F% < 1 OR F% > 6 


BEEP 

LOCATE 12, 26: PRINT "Enter desired frequency (1-6)"; : F% = 
VAL (INPUTS (1) ) 
WEND 
SELECT CASE F% 

CASE F% = 1: RclPanel$ = "0F010PFF" "PANEL 1 40 

CASE F% = 2: RclPanel$ = "OFO20DFF" "PANEL 2 100 

CASE F% = 3: RclPanel$ = "OFO30DFF" "PANEL 3 315 

CASE F% = 4: RclPanelS = "OFO40DFF" "PANEL 4 1000 

CASE F% = 5: RclPanel$ = "0F050PFF" "PANEL 5 1500 

CASE ELSE: RclPanel$ = "0F060DFF" 'PANEL 6 All of them 
END SELECT 
COLOR 30, 0: LOCATE 16, 16 'Blinking yellow on black 
PRINT "Performing initial setup."; : COLOR 14, 0 


PRINT "This will take a few seconds." 
CALL IBSAD(SD375%, 115) 
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CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 


IBWRT (SD375%, 
IBSAD(SD375%, 
IBWRT(SD375%, 
IBWRT(SD375%, 
IBWRT(SD375%, 
IBWRT(SD375%, 
IBSAD(SD375%, 
IBWRT (SD375%, 
IBSAD(SD375%, 
IBWRT (SD375%, 


"MPB(A)/" + CHR$(13)) 
111) 

RclPanel$) 
"1416010000000BFF") 
"1516010000000BFF") 
"1F1306000BFF") 

115) 

"MRI(A)/" + CHRS(13)) 
107) 

"330DFF") 


'MAKE PUSH BUTTONS ASCII 
'PUSH BUTTON GROUP 15 
'PANEL $_ RCL 

'CH A 1000mv/EU 

'CH B 1000mv/EU 

'dB, AVG T = 60 SEC 
"GROUP 19 CONTROL 

"MAKE RANGE INPUTS ASCII 
"GROUP 11 RANGE INPUTS 

0 = 2 volt, 0- 2 Khz 


"kkkkkkkkkkkkkkkk ready to start data acquisition *xxxxkkkkkkkkkkkkkkkkk 


END SUB SUB GetPrecisionData STATIC 

SHARED SD375*, PrecisionData AS STRING * 6402, DataDateS, T1! 
SHARED ChA AS RecordTypel, ChB AS RecordTypel, SampleTime$ 
SHARED NumberOfRecordsi, NumberOfRecords2, NumberOfRecords3 


Cell(1) = 16 ' 40 Hz ----- vade 
Cell(2) = 40 Dee o 
Cell(3) = 126 ' 315 Hz ----- _---===) CHANNEL A 
Cell(4) = 200 "1000 Hz ----- cn 
Cell(5) = 300 '1500 Hz ----- 
Cell(6) = 416 ' 40 Hz ----- o 
Cell(7) = 440 ' 100 Hz ----- mr 
Cell(8) = 526 5 :5—-.-. 2 ---===) CHANNEL B 
Ce11(9) - 600 10005 TT . --- 
Се11(10) = 700 '1500 Hz ----- 
LOCATE 1, 1: PRINT "Getting Precision data..."; SPACES(45) 
SampleTime$ = TIMES 
TIMER OFF 


CALL IBSAD(SD375%, 100) 


CALL IBRD(SD375%, PrecisionData) ' 3.2 sec 
TIMER ON: T1! = 0 
LOCATE 1, 1: PRINT SPACES(26); "...data received..."; SPACES(15) 


CALL IBSAD(SD375%, 115) 
CALL IBWRT(SD375%, "MPB(A)/" + CHR$(13)) 
CALL IBSAD(SD375%, 111) 


"CONTROL GROUP 
'MAKE PUSH BUTTONS ASCII 
'PUSH BUTTON GROUP 15 


CALL IBWRT(SD375%, "31474846FF") 
'starts the cycle timer here 


MI, STOP, ERASE, START 


FOR 1 = 1 TO 10 
First% = 8 * (Cell(I) - 1) + 3 
cell (I) 
CellStringS = MIDS(PrecisionData, First, 8) 
"Compute BYTES 1 and 2 
FOR J = 1 TO 2 
K=2* J-1 
H1$ - MIDŞ(CellString$, K, 1) 
H2$ = MID$(CellString$, K + 1, 1) 
IF HIŞ < "A" THEN D1% = ASC(H1$) - 48 ELSE D1% = ASC(H1$) - 55 


'beginning of bytes for 
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IF H2$ < "A" THEN D2% = ASC(H2$) - 48 ELSE D2% = ASC(H2$) - 55 
Byte%(J) = D1% * 16 + D2% 

NEXT J 

‘BYTE 3 is always zero 

'Compute BYTE 4 
H1$ = MID$(CellString$, 7, 1) 
H2$ = MID$(CellString$, 8, 1) 
IF H1$ € "A" THEN D1% = ASC(H1$) - 48 ELSE D1% 
IF H2$ < "A" THEN D2% = ASC(H2$) - 48 ELSE D2% 
Byte%(4) = D1% * 16 + D2% 
Real! (I)=MakeReal! (Byte%(1) ,Byte%(2) ,Byte%(4))/451'for hanning window 
LOCATE 13 + I, 1 
Log10! (I) = LOG(Real!(I)) / LOG(10#) 


ASC (HIS) TEES 
ASC (H25) 02055 


PRINT I; TAB(10); Real! (I); TAB(26); " Volts "> : COLOR 9090 
PRINT 10 * Log10!(I); TAB(48); "dB ": COLOR 14, O 
NEXT I 
LOCATE 1, 1: PRINT SPACES(3); "...writing to files."; SPACES(30) 
"put first 5 values in channel A file, second 5 in Channel B file 
ChA.Freql = Real! (1): ChB.Freql = Real! (6) 
ChA.Freq2 = Real! (2): ChB.Freq2 = Real! (7) 
ChA.Freq3 = Real! (3): ChB.Freq3 = Real! (8) 
ChA.Freq4 = Real! (4): ChB.Freq4 = Real! (9) 
ChA.Freq5 = Real! (5): ChB.Freq5 = Real! (10) 
ChA.TimeO = SampleTimeS: ChB.Time0 = SampleTime$ 
ChA.DataDate = DataDateS: ChB.DataDate = DataDateS 


"write to file 

RecordNumber2 = RecordNumber2 + 1 

RecordNumber3 = RecordNumber3 + 1 

PUT $2, RecordNumber2, ChA 

PUT #3, RecordNumber3, ChB 

WRITE #5, MIDS (PrecisionData,3,3200) ‘seq file for raw ch A data 
WRITE $6, MIDS(PrecisionData, 3203, 3200) 'raw data for Ch B 
LOCATE 1, 1: PRINT "...writing completed. 2 
END SUB 


SUB AskIfPrintoutWanted 
CLS : LOCATE 12, 20 
INPUT "Do you want a printout of the data files (y/n)? ", AS 
DO WHILE UCASES(AS) <> "Y" AND UCASES(AS) <> "N" 
LOCATE 12, 20 
INPUT "Do you want a printout of the data files (y/n)? ", AS 
LOOP 
CLOSE #1: CLOSE #2: CLOSE #3 
IF UCASES(AS) = "Y" THEN CALL PrintFiles 
END SUB 


DEFINT A-Z 

SUB PrintFiles 

SHARED ChA AS RecordTypel, ChB AS RecordTypel 
SHARED Spectrum AS RecordType2, Volt(), Level() 
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OPEN "ChanA.dat" FOR RANDOM ACCESS READ AS #2 LEN LEN (ChA) 

OPEN "ChanB.dat" FOR RANDOM ACCESS READ AS #3 LEN LEN (ChB) 

OPEN "Voltage.DAT" FOR RANDOM ACCESS READ AS #1 LEN = LEN(ChA) 
NumberO0fRecords1=LOF (1) /LEN(Spectrum) 'Determines $ records in each file 


NumberOfRecords2=LOF (2) /LEN (ChA) 'so that existing files will not 
NumberOfRecords3=LOF (3) /LEN (ChB) 'be written over later. 

LOCATE 7,1:PRINT "Number Of Records in Spectra.DAT is ";NumberOfRecords1 
TS =" DATE VCR TIME 40 Hz 100 Hz 315 Hz 1000 Hz 1500 
Hz" 

L$="' I A E ———. A A = —.-..—.—...— ......... —.... 


MM-DD-YYYY #-A HH:MM:SS #. 888888 +. HHHH #. HHHH А АНАНЫН 
E. EHHH 


' k x x Print Spectra.DAT file K k K X k k k X k k k k k Kk k k Kk k k k k k k k k k k k k k k k k k k x x K 


PRINT #4,"Spectra.DAT data file ";"current time is ";TIMES; " ";DATES 
FOR RecordNumber = 1 TO RecordNumberl 

GET $1, RecordNumber, Spectrum 

Volt (RecordNumber) = Spectrum.Voltage 

Level (RecordNumber) = Spectrum.Level 
NEXT RecordNumber 


US = "të fit GR. HÅ Hti# İ.İH 3157 BHH 88" 
FOR I = 1 TO 50 ‘print to screeen until tested 
J = I + 50 


PRINT USING U$; 5 * I; Volt(I); Level(I); 5 * (J); Volt(J); Level(J) 
NEXT I 
PRINT CHRS (12) 
FOR I = 101 TO 150 
J = I + 50 
PRINT USING U$; 5 * I; Volt(I); Level(I); 5 * (J); Volt(J); Level(J) 
NEXT I 
FOR I = 201 TO 250 
J =I + 50 
PRINT USING U$; 5 * I; Volt(I); Level(I); 5 * (J); Volt(J); Level(J) 
NEXT I 
PRINT CHRS (12) 
FOR 1 = 301 TO 350 
J = I+ 50 
PRINT USING US; 5 * I; Volt(I); Level(I); 5 * (9g); Volt(J): Level(J) 
NEXT I 
PRINT CHRS (12) 
Goto bypass: 
t k k k k k k à k k x K Print Channel A data file K k x x k k k k k k k k k k k k k k k k k k k k k k K 
PRINT $4, "Channel A File data": PRINT #4, T$: PRINT $4, LS 
FOR RecordNumber = 1 TO NumberOfRecords2 
GET 2, RecordNumber, ChA 
PRINT #4, ChA.DataDate; " "; ChA.ChannelID; " "; ChA.Time0; " "; 
PRINT #4, USING " #.######°° °°"; ChA.Freql; ChA.Freq2; ChA.Freq3; 
ChA.Freq4'; ChA.Freq5 
NEXT RecordNumber 
PRINT $4, CHRS(12) 'form feed 


' k k x x x k x x x k x x x x Print Channel B data file K K K K k k k k k k k K k k k k k k k k XK k k k k k k k k k 


PRINT #4, "Channel B File data": PRINT $4, TS: PRINT #4, LS 


gi 


FOR RecordNumber = 1 TO NumberOfRecords3 

GET $3, RecordNumber, ChB 

PRINT #4, ChB.DataDate; ChB.ChannelID; " "; ChB.Time0; " "; 

PRINT #4, USING " #.######°°°°"; ChB.Freql; ChB.Freq2; ChA.Freq3; 
ChB.Freq4'; ChB.Freq5 
NEXT RecordNumber 

Bypass: 
END SUB 


SUB CheckStatus 

SHARED SD375%, Averaging%() 

CALL IBSAD(SD375%, 115) 

CALL IBWRT(SD375%, "MST(B)/" + CHRS(13)) 


CALL IBSAD(SD375%, 111) 'SD 375 must be in BINARY Byte coding 
CALL IBRDI(SD375%, Averaging%(), 2) '"1" if avg in progress (p.4-23) 
TO! = TIMER: Delay! = 0! 
DO UNTIL Delay! > .4 'Waits .4 SEC before next status check 
Delay! - TIMER - TO! 
LOOP 
END SUB 


SUB ConvertRawData 
SHARED Spectrum AS RecordType2, RawData AS STRING * 3200 
SHARED Filel$, Sample% 

'*** check to see if random access file already exists **** 

CLS: PRINT "Already have random file to break up (Y/N)?": AS - INPUTS (1) 
IF UCASES(AS) = "Y" GOTO StopHere '*** skip file creation if 'yes' *** 
'** get name of ASCII input file (created by SD 375 transfer process) 
INPUT "Complete name of input file (FILENAME.EXT)?", InputFileS 
PRINT "If input file does not exist this program will hang up." 
LOCATE 12, 20 
PRINT "Enter [A] to abort or [C] to continue.": AS = INPUTS(1) 
'*** test for correct response Xxx 
DO WHILE UCASES(AS) © "A" AND UCASES(AS) C) "C" 

LOCATE 12, 20: BEEP: BEEP 

PRINT "Enter [A] to abort or [C] to continue.": AS = INPUTS (1) 
LOOP 
IF UCASES(AS) = "A" THEN GOTO StopHere: '* abort program if "A" entered 
'* convert ASCII file to random access file * 
CLS : LOCATE 12, 20: PRINT "Converting data..." 
OPEN InputFileS FOR INPUT AS #10 
OPEN "SPECTRA.DAT" FOR RANDOM ACCESS WRITE AS #17 LEN = LEN(Spectrum) 
'* SPECTRA.DAT is final product (used to create additional file later) * 
RecordNumber = 0 
DO UNTIL EOF(10) 

INPUT #10, RawData '* read 1 string of ASCII HEX data (400 groups of 8) 
FOR I = 1 TO 400 "250 stops at 1250 Hz (there are 400 total for 2KHz) 
First% = 8 * (I - 1) + 1 ‘beginning of bytes for cell I 
CellStringS = MIDS(RawData, First%, 8) 

'** Compute BYTES 1 and 2 ** 
FOR J = 1 TO 2 
K= 2* J - 1 
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H1S = MIDS(CellStringS, K, 1) 
H2$ = MIDS(Cellstring$, K + 1, 1) 
IF H1$ < "A" THEN D1% = ASC(H1$) - 48 ELSE D1% 
IF H2$ < "A" THEN D2% = ASC(H2$) - 48 ELSE D2% 
Byte%(J) = D1% * 16 + D2% 
NEXT J 
'** BYTE 3 is always zero ** 
'** Compute BYTE 4 ** 
H1$ = MIDS(CellString$, 7, 1) 
H2$ = MID$(CellString$, 8, 1) 
IF H1S < "A" THEN D1% ASC(H1$) - 48 ELSE D1% 
IF H29 < "A" THEN D2% ASC(H2$) - 48 ELSE D2% 
Byte%(4) = D1% * 16 + D2% 
FSV=1 'Full Scale Value for SD375 (=1 7/27/87 data,=2 7/28/87 data) 
Spectrum.Voltage = FSV * (MakeReal!(Byte%(1),Byte%(2),Byte%(4))/451) 
Spectrum.Freq = 5 * I ‘volts/451 removes gain from Hanning window 
Spectrum.Level = 10 * LOG(Spectrum.Voltage) / LOG(10#) ' MAG. not 
POWER 
RecordNumber = RecordNumber + 1 
PUT #17, RecordNumber, Spectrum 
NEXT I 
LOOP 
CLOSE #15: CLOSE #16: CLOSE #17 '*** close all open files *** 
PRINT "Do you want to continue? (Y/N) ": AS = INPUTS(1) 
IF UCASES(AS) = "N" THEN GOTO EndProg 


ASC(H19) — 55 
ASC(H2S) $9455 


ASC (H1S) 55 
ASC(H2$) - 55 


StopHere: 
'*** get source filename (created above) for creation of small files *** 
INPUT "Enter name of random Input file (file MUST exist): ", InputFileS 


OPEN InputFile$ FOR RANDOM ACCESS READ AS #7 LEN = LEN(Spectrum) 
NumberOfSamples = (LOF(7) MV LEN(Spectrum)) 1 400 ‘'250 is for 1250 Hz 
MakeChoice: 

CLS : COLOR 14, O: LOCATE 10, 15 

PRINT "Number of samples in "; InputFile$; " is "; NumberOfSamples 
LOCATE 12, 15 ' *** separate file for each sample *** 

PRINT "If you want to make one-sample spectra files, press [A]" 
LOCATE 13, 15 ' *** separate file for each frequency *** 

PRINT "If you want to make single-frequency files, press [B]" 
LOCATE 16, 15 

PRINT "Enter your choice [A], [B], or [E] to exit.": AS - 
UCASES ( INPUTS (1)) 


COLOR 9, O 
SELECT CASE AS 
CASE "A" '"X***** CREATE Å SEPARATE FILE FOR EACH SAMPLE ****** 
Sample% = 1 
DO 


CLS : LOCATE 1, 15: 

INPUT "Enter full name of output file for magnitude: ", Filel$ 
LOCATE 14, 15 : INPUT "Enter sample number:", Sample% 

First = 400 * (Sample% - 1) + 1 

OPEN FileiS FOR OUTPUT AS #8 

VariableNameS = LEFT$(File1$, 3) 


73 


PRINT $8, VariableName$ 
FOR I = First TO First + 399 '249 FOR 1250 HZ 
GET #7, I, Spectrum 
PRINT #8, Spectrum.Voltage 
NEXT I 
CLOSE #8 
LOCATE 16, 15 
PRINT "Do you want to make another file? (Y/N) ": Q$ = INPUT$ (1) 
Sample% = Sample% + 1 
LOOP UNTIL Sample% > NumberOfSamples UCASES (QS) = "N" 
CASE "В" ' *xxxxx CREATE A SEPARATE FILE FOR EACH FREQUENCY **xxxx 
DO 
CLS : LOCATE 10, 1: 
INPUT "Enter name of output file for magnitude:", FS 
LOCATE 14, 10: COLOR 9, 0 
DO 
INPUT "Enter frequency (MULTIPLE OF 5, 5 =< F <= 2000): ", F% 
CLS : LOCATE 14, 10: COLOR 14, O 
PRINT "Frequency entered is "; F%; " Hz. Is this correct? (Y/N) " 
AS = INPUTS (1) 
IF UCASES (AS) = "N" THEN | 
INPUT "Enter frequency (MULTIPLE OF 5, 5 =< F <= 2000):", F% 
END IF 
LOOP UNTIL UCASES(AS) = "Y" 
COLOR 9, O "blue on black 
S% = 0 "keeps track of sample number 
Bin% = F% ÀN 5 'bin number of each sample 
OPEN FS FOR OUTPUT AS #8 
FOR I = 1 TO NumberOfSamples 
Record’ = (I - 1) * 400 + Bins 
GET #7, Records’, Spectrum 
PRINT #8, Spectrum.Voltage 
NEXT 1 
CLOSE $8 
LOCATE 16, 15: COLOR 14, 4: 'yellow on red 
PRINT "Do you want to make another file? (Y/N) ": Q$ = INPUTS (1) 
COLOR 14, 0 "yellow on black 
PRINT "Number of values in file is "; NumberOfSamples 
LOOP UNTIL UCASES(QS) = "N" 
CASE "E": GOTO EndProg 
CASE ELSE: BEEP: BEEP: BEEP: GOTO MakeChoice 
END SELECT 
EndProg: 
END SUB 


SUB ResetClock 

TIMER OFF 

CLS : COLOR 9, 4 

PRINT "You must enter the current time for the computer system clock." 
COLOR 14, O 

DO: INPUT "The Hour is (00 - 24) ", Hour$ 

LOOP UNTIL VAL(HourS) >= 0 AND VAL(HourS) <= 24 
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PRINT 

DO: INPUT "The minutes are (00 - 59)", Minutes$ 

LOOP UNTIL VAL(Minutes$) >= 0 AND VAL(Minutes$) <= 59 
PRINT 

DO: INPUT "The seconds are (00 - 59) ", SecondsS$ 

LOOP UNTIL VAL(Seconds$) >= O AND VAL(Seconds$) <= 59 
PRINT 


TIMES = Hour$ + ":" + Minutes$ + ":" + Seconds$ 
END SUB 
DEFSNG A-Z 


FUNCTION MakeReal! (B1%, B2%, B4%) 
MakeReal! = (256 * B1% + B4% - 8192) * (2 * (B2% - 128)) / 8192 


END FUNCTION 


T3 


APPENDIX B 


AMBIENT NOISE LEVEL PLØTS 
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Figure C2. Visual satellite photograph taken at 12 


55 on 25 July 87. 
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